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Abstract 

— s ■ 

We study the equilibrium spin configuration of the 2D Hubbard model at 
low doping, when a long-range magnetic order is still present. We use the 

^ ■ 

' spin-density-wave formalism and identify three different low-doping regimes 

o' 

depending on the value of z = 4Ux2D where \2D is the Pauli susceptibility of 
holes. When z < 1, the collinear antiferromagnetic state remains stable upon 



low doping. As candidates for the ground state for z > 1 we first examine 

a" 

the planar spiral phases with the pitch Q either in one or in both spatial 



directions. Mean-field calculations favor the spiral (tt, Q) phase for 1 < z < 2, 
and (Q, Q) phase for z > 2. Analysis of the bosonic modes of the spiral state 



^ ■ shows that the (Q, Q) state has a negative longitudinal stiffness and is un- 



stable towards domain wall formation. For the (it, Q) state, the longitudinal 
stiffness is positive, but to the lowest order in the hole concentration, there is 
a degeneracy between this state and a whole set of noncoplanar states. These 
noncoplanar states are characterized by two order parameters, one associated 
with a spiral, and the other with a commensurate antiferromagnetic ordering 
in the direction perpendicular to the plane of a spiral. We show that in the 
next order in the hole concentration this degeneracy is lifted, favoring non- 
coplanar states over the spiral. The equilibrium, noncoplanar configuration 
is found to be close to the Neel state with a small spiral component whose 



amplitude is proportional to the square root of the hole concentration. These 
findings lead to a novel scenario of spin reorientation upon doping in Hubbard 
antifer r omagnets . 
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I. INTRODUCTION 



Magnetic properties of the Cu0 2 layers in high temperature superconductors have been 
recently attracting an intense interest as magnetism is possibly a major contributor to the 
mechanism of superconductivity [II] . There are numerous reasons to believe that the mag- 
netic properties of weakly doped cuprates are quantitatively captured by the effective 2D 
theory for one degree of freedom per CuOi unit which is provided by a one-band Hubbard 
model 

H = -Y, kj a ta a j,a + n^n iA (1) 

Here a is a spin index, n = a^a, and ti j is a hopping integral which acts mainly between 
nearest (t) and next-nearest (f) neighbors. We will assume that t' is negative [0. 

At half-filling, the ground state of the 2D Hubbard model exhibits a long range commen- 
surate Neel order provided t' is not very large. It has been known for many years H,|J that 
holes introduced into a commensurate antiferromagnet give rise to a long-range dipolar dis- 
tortion of the staggered magnetization. In 2D this effect was studied in detail by Shraiman 
and Siggia |J. They found that in the simplest scenario, the dipolar distortion leads to a 
spiral spin configuration with the momentum (n, Q) at any nonzero doping. The incommen- 
surate (7r, Q) phase was also obtained in the early perturbative studies of the Hubbard model 
with small U ||, and in several other mean- field [jTj-fl and self-consistent ||10|| calculations. 



In this paper, we use the spin- density- wave (SDW) approach and study the structure 
of magnetic correlations in the Hubbard model at small but finite doping when long-range 
magnetic order is still present. We will show that depending on the strength of the interaction 
between the holes, three different solutions of the Hubbard model at low doping are possible: 

(i) commensurate Neel phase at small interaction (precise criteria will be derived below), 

(ii) phase separation at sufficiently strong interaction, and (iii) intermediate homogeneous 
incommensurate phase which, however, differs from the planar spiral suggested by the mean- 
field analysis. This new incommensurate phase is noncoplanar in the spin space, and its 



3 



properties are much closer to the properties of a commensurate (it, tt) spin configuration 
than those of a planar spiral. 

The analysis presented here is related to other works on incommensurate magnetic phases 
at finite doping. Previous mean- field studies of the Hubbard and t — J models PJ7|-p|JTT| 
have focused on the three configurations with ordering momenta (vr,7r), (tt, Q) and (Q,Q), 
and have shown that in a certain range of parameters, any of these configurations can have 
energy lower than that of the other two states. Our energy analysis is consistent with their 
results. These mean-field studies also found that for short-range repulsion, the (Q, Q) phase 
is likely to be unstable towards domain walls formation ||||, while the (tt, Q) phase is 
stable. Shraiman and Siggia |12| developed a macroscopic theory of the bosonic excitations 
in the (n, Q) spiral phase. Surprisingly, to the lowest order in density they found a peculiar 
degeneracy in the ground state energy for the planar spiral state and for a whole set of 
noncoplanar magnetic configurations with the plane of the spiral varying in space. This 
degeneracy is also present in our microscopic calculations reported below. However, the 
further assumption of Shraiman and Siggia that the next-order terms in doping concentration 
stabilize the planar spiral state is inconsistent with our microscopic calculations which favor 
noncoplanar spin configuration for the same range of parameters as in JT2I . 



The paper is organized as follows. In Sec [H] we consider the mean-field theory of the 
spiral phases in the Hubbard model and compare the ground state energies of different 
phases. In Sec [Hj we compute the dispersion of the bosonic excitations in the (tt, Q) and 
(Q, Q) spiral phases to the lowest order in hole density and find that the (Q, Q) spiral has 
negative longitudinal susceptibility, while the (tt, Q) spiral has an infinite number of zero 
modes. We will identify the set of magnetic states which are degenerate in energy with 
the planar (tt, Q) spiral. In the next section, we show that the degeneracy is lifted by the 
next-to-leading order terms in hole density. We compute the ground state energy and find 
the equilibrium spin configuration at a finite doping. We then discuss the properties of this 
equilibrium state, in particular, the behavior of the dynamical spin susceptibility. Finally, 
Sec IVTI states our conclusions. 



II. THE SPIRAL PHASES 



In this section, we describe spin-density wave (SDW) calculations for the Hubbard model 
at and near half-filling. Let us start with the commensurate Neel state. The SDW approach 



the results of these studies. At half-filling, the fermionic spectrum consists of conduction 
and valence bands separated by the energy gap A = U(S Z ). The dispersion relation for the 



valence fermions is Ef. = —E~ + e + where E- = J A 2 + (e~) 2 and = — 2t(co$k x + cos k y ), 

e + = A\t' cos k x cos k y \. This dispersion has a maximum at four points (±7r/2, ±7r/2) in the 
center of each of the edges of the magnetic Brillouin zone, provided that t' is not too large. 
In the neighborhood of these points, E% can be presented as Ef = — A + p%/2m\\ + p]_/2m±. 
Thus near (7r/2,7r/2), we have p\\ = (k x — k y )/2, p±_ = {k x + k y )/2), m\\ = (4t') _1 , m± = 
(4J - 4|t'|)- 1 where J = At 2 /U is the inverse bandwidth. For both La— and Y— based 
materials, t' ~ J and therefore both masses scale as 1/J. Note, that the minimum at 
(7r/2,7r/2) is rather robust — even if t' = and the mean-field spectrum is degenerate 
along the whole edge of the magnetic Brillouin zone, the actual dispersion still has four 
minima at (±7r/2, ±7r/2) due to quantum fluctuations (not included in our present SDW 
treatment) 0,[18|. 

At finite doping, the chemical potential moves into a valence band, and the states near 
the maximum of Ef become empty. These states are often referred to as hole pockets. For 
a commensurate spin ordering, all four pockets become equally occupied. Performing then 
simple calculations, we obtain that in the presence of holes, the ground state energy of the 
(-7T, 7r) phase changes to 

= nx " (2 ) 

We will also need the expression for the static magnetic susceptibility of the (tt, it) state. 
In the SDW formalism, the susceptibility is obtained by summing up the series of bubble 
diagrams, and the result is 



for this state has been discussed several times in the literature |7,13-16| and we will use 
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x + (<?) 



where for A > t, t' 
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(4) 



Prime at the summation signs indicate that the summation is over the magnetic Brillouin 
zone. Near Q = (tt,tt), the static transverse susceptibility should obey the a hydrodynamic 
relation Jl9| XstT(<?) = 2iVg/(p s (g — Q) 2 ), where N is the on-site magnetization, and p s is 
the spin-stiffness. At half-filling only bubbles containing valence and conduction fermions 
are allowed, and performing simple calculations, we obtain p s = p° = J(l — 2(t'/t) 2 )/4, 
where J = 4t 2 /U. At finite doping, there is also a contribution to the stiffness from the last 
term in (|j) which contains only valence fermions. This last contribution is proportional to 
the Pauli susceptibility, which in two dimensions does not depend on carrier concentration. 
As a result, in the SDW approximation, the spin stiffness acquires a finite correction at an 
arbitrarily small deviation from half-filling P,p|,|I6[| 



Ps = - z), 



(5) 



where 



z = 4Ux 



P ,,,ii _ 2u} /fnjm\\ 



2D 



TV 



(6) 



We see that if z < 1, the Neel state remains stable at finite doping, while if z < 1, it 
becomes unstable, and we have to consider incommensurate states as possible candidates 
to the ground state. In the naive mean-field calculations discussed so far, we have z ~ 
U/J ^> 1, which implies that commensurate state becomes unstable immediately upon 
doping. However, we have explicitly verified that all our conclusions are valid for arbitrary 
ratios of t/U and t'/t provided that magnetic order at half-filling is commensurate, and 



the hole pockets are located at (±7r/2, ±7r/2). We therefore will consider z as an imput 
parameter which can, in principle, have any value. Some rigorous results about the dispersion 
of bosonic excitations at arbitrary t,t', U will be presented in appendix B. 

It is worth emphasizing that even in the large U limit, the value of z can in fact be 
of the order of unity. The point is that the mean-field results at large U must be taken 
with caution in view of strong self-energy and vertex corrections which both contribute 
powers of U/tS. Self-consistent consideration of these corrections indicate that they do not 
change the momentum dependence of the vertices at small q — Q, but reduce the overall 
scale of the effective interaction between holes from U to U e ff which is of the order of the 
bandwidth J P JIEkpOf . This in turn implies that z is in fact simply a number, independent of 
U/t. Furthermore, striktly speaking, instead of U e s we have to consider the total scattering 
amplitude of two holes, T. In two-dimensions, 

T = ^ (7) 

Hence T, and consequently, z vanishes logarithmically as the hole concentration, 5 oc y/p F , 
tends to zero. This always makes the collinear antiferromagnetic state stable at very low 
doping [IB]. However, the range of doping where logarithmic corrections are important is 



likely to be very small, and in this paper we simply set T = U e s and consider z = 0(1) as a 
doping-independent parameter. For simplicity, throughout the paper, we will focus on large 
U and present the results for the quantities, not related to Pauli susceptibility, only to the 
leading order in t/U. 

We now consider incommensurate spin configurations at finite doping. Let us first fo- 
cus on the two simplest candidates: the spiral states with the ordering vectors (tt, Q) and 
(Q,Q) HJ9|-pl]|. As before, we will label the ordering momentum as Q (Q can be either 
(tt, Q) or (Q, Q)). The spin order parameter now has two components, and in terms of the 
fermionic operators is expressed as 

k 



SZq= 5>L«*-«,t> = ( S qT ■ (8) 
fe 

For definiteness, we choose the order parameter to be in the XY plane. Without a loss of 
generality one can also choose Sq to be real. The real space spin configuration described by 
(§) is then = Sq cos(QR), and = Sq sin(QR), where R denotes the lattice site. 

The SDW calculations proceed in the same way as before: one has to decouple the 
interaction term in ([TJ) using (^) and diagonalize the quadratic form. Performing the com- 
putations, we obtain 

Ef = t + ±E_ 1 (9) 



where as before EL = y A 2 + (e_) 2 , but now 

_ e k+Q/2 + e k-Q/2 _ e k+Q/2 ~ e k-Q/2 n m 

e+ ~ 2 ,e -~ 2 • ( U) 

Note that for XY-ordering, these new electron states appear as hybridization of the original 
electrons with opposite spins and thus contain no spin labels. In this situation, there is no 
doubling of the unit cell, and the summation over momenta in (]TT[) is extended to the whole 
first Brillouin zone (— tt/cl < k x>y < ir/a). 

The gap A is related to the parameters of the Hubbard model via the self-consistency 
condition 

1 W , 1 in, 

V N jf 2,/A* + El 



where the summation goes over the momenta of occupied states. 

Consider now specifically the spiral (tt, Q) state. The mean-field fermionic spectrum for 
this state is not symmetric with respect to the reflection k y — > —k y , although it is still 
symmetric with respect to k x —>■ —k x . Consequently, the minima at (±7r/2,7r/2) have lower 
hole energy than those at (±7r/2, —tt/2). 

Suppose that the concentration of holes filling the lower energy pocket is X\, and the 
higher energy pocket — x<i [x\ + £2 — x). Simple calculations then show that the ground 
state energy of the (tt, Q) phase is given by 
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E^ = -tq( Xl -X2) + i Sr + , (12) 

4A 2 y /m±m\\ 

where q = q y = n — Q. The first term reflects the decrease in total energy due to the e + term 
in the spectrum, the second term results from the redistribution of the energy levels below 
the Fermi level, and the last term corresponds to the increase of the total energy due to 
unequal occupation of pockets. Minimizing the total energy of the (it, Q) state with respect 
to q we obtain to the lowest order in hole concentration 

q = j(xi - x 2 ). (13) 

For the energy of the (it, Q) state, we then have 

= ^-^ (i - z)+ ™ 2 . (14) 

We see that the ground state energy of the (ir, Q) phase becomes smaller than that of the 
(it, 7r) phase at z > 1, exactly when the stiffness for the (if, it) phase becomes negative. In 
the latter case, we also have X\ — x and x 2 = 0, which implies that only two out of four 
pockets are occupied, and q = (U/t)x. A similar analysis was performed in |J. 

Finally, consider the (Q, Q) phase. Now the fermionic spectrum is not symmetric with 
respect to either the k x — > —k x or k y — > —k y reflections. The point (tt/2,tt/2) becomes the 
only absolute minimum of the hole spectrum. Suppose that the concentration of holes filling 
the lowest energy pocket is X\, the two intermediate energy pockets - x%, and the highest 
energy pocket - x 3 (x\ + 2x 2 + x 3 = x). Then the ground state energy of the (Q, Q) phase 
is given by 



E (q,q) = _ 2tq{xi _ X3) + t ^ + ^M ±^^ L _ 

The inverse pitch of the spiral q = it — Q is found by minimizing the energy: 



(15) 



q=j(xi-x 3 ). (16) 



The total energy then assumes the following form: 



E (Q,Q) = ( t Xi _ x \2(\ - z ) + ( x - 2x 2 ) 2 + 4x1) . (17) 

It is immediately obvious, that if z < 1, then the lowest possible energy is achieved when 
x\ = X2 = X3, and q = 0, i.e. in the (tt, tt) phase. If z > 1, then 23 = 0. Minimizing the 
energy with respect to X\, we find 

4 

xi = x (18) 

6 — 2 

when z < 2, and 

xi = x (19) 

when 2: > 2. In the former case only the highest energy pocket has no holes, and the total 
energy is equal to 

E (Q,Q) = 2ttx 2 2-£^ 
y / mlmj[6 — 2; 

It is straightforward to see that it is always higher than the energy of the (tt, Q) phase, given 



by eq. (|H ) . In case of z > 2, only the pocket with the lowest hole energy is occupied, and 
the total energy is 



TlX 2 



E (Q,Q) _ — _ ( 2 - z l (21) 

Comparing (0) ( P^|) and (PT|), we observe that the Neel state is the minimum for z < 1, the 
(7r, Q) spiral phase has the lowest energy at 1 < z < 2, while the (Q, Q) state has the lowest 
energy at z > 2. 

Observe however, that in the latter case, d 2 E^'^ /dx 2 ~ (2 — z) is negative. On general 
grounds, this result suggests that the homogeneous solution is unstable. We will see in the 
next section that the longitudinal stiffness for the (Q, Q) state is in fact negative - this 
will be another argument in favor of an inhomogeneous ground state. On the other hand, 
for the (tc,Q) phase in its region of stability (1 < z < 2) we have, dE^'^/dx 2 > 0, i.e., 
homogeneous solution is stable. 
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III. COLLECTIVE EXCITATIONS 



Now we proceed to examining the stability of the spiral states by considering collective 
bosonic excitations. It follows from general considerations that in a spiral phase one should 
have a Goldstone mode with one velocity related to the spin rotation in the plane of the 
spiral, and two Goldstone modes with another velocity related to the rotations of the plane of 
the spiral around X and Y axes [BTJ. The former Goldstone mode results in the divergence 
of the total static susceptibility x + ~(l) & t the wave vector q = —Q, while the latter lead to 
divergences of the total x zz {l) a t the wave vectors ±Q. 

The spectrum of collective excitations is determined by the poles of dynamic suscepti- 
bility 

x£ = i/ dte^(TSl(t)Si q (0)) (22) 

Here S l represents either one of the three spin densities, S + , S~ , S z , or the charge density 
p. For the (n, ir) state, fluctuations in transverse spin channels are completely decoupled 
from fluctuations in the density and longitudinal spin channels. Dynamical susceptibility 
is then a 2 x 2 problem. For planar spiral states, however, all four channels are coupled, 
and the dynamical susceptibility has to be found by solving a set of four coupled Dyson 
equations. Doing the standard SDW manipulations, we obtain that the poles of dynamical 
susceptibility are given by solving D(q,u) = where D(q,u) is the determinant of a 4 x 4 
matrix given by 



~ U X q ^( q+2 Q) 1 ~ U X g +2Q-(q+2Q) ~ U ^X^Q -(q+2Q) U ^X^.Q-^ q+2 Q) 
-Uy/2Xq*_( q+ Q) ~U V^Xg+ 2 Q ,-{q+Q) 1 _ 2U K+Q ,-(q+Q) 2UX q+Q-(q+Q) 
V ~ Uy//2 ~ X t-(Q+Q) ~ Uy//2 ~ X (q+2Q),-(q+Q) - 2U X q P + Q^ q+ Q) 1 + 2U X q P + Q _(,+Q) J 

The expressions for the irreducible susceptibilities are presented in the appendix A. A similar 
expression for D(q,u) was recently obtained by Cote and Tremblay |22[] in their SDW anal- 
ysis of collective excitations in 2D Hubbard model on a triangular lattice at half-filling For 
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our present consideration, it is essential that at zero frequency, x q -( q +Q) = Xg+ 2 Q -{q+Q) = 
Xq+g -(q+Q) = at any q, and therefore static transverse spin fluctuations decouple from the 
longitudinal spin fluctuations and density fluctuations (we remind that spin ordering is in 
the XY plane). We now consider the transverse and longitudinal fluctuations separately. 



A. Longitudinal spin fluctuations 

Consider first the solution for the density and longitudinal spin fluctuations. At q = —Q, 
the evaluation of the expressions in the appendix A yields at zero frequency xVq = Xq P — 
X^q = z/8U, x^J,_g = = (z- 4)/8U, and X % Q - ^ = E fe l/(2E k ) ee l/U. 
Elementary manipulations then show that there is indeed a Goldstone mode at q — —Q. 
Expanding around this point, we obtain after straightforward but lengthy calculations, that 
to quadratic order in q + Q 

/ ~\2 ( ' Zui\ 



D{^) = -{q + Qf[l-^--. (23) 

The q + Q term comes from the expansion of 1 — UXqZ q + UXg*r q+2 Qy an d ^ term comes 
from x +z an d X~ z ■ We also introduced z w = 4Ux2d( uj ) where X2d(^) is the susceptibility of 
a 2D Fermi gas at finite frequency 



X2D{U) = 



1/2- 



2tt 

where to the leading order in the hole density, 



1 - 



S 2 



(24) 



and pf ~ y/x is the Fermi momentum of holes. At zero frequency we indeed have = z. 
The explicit expression for the total longitudinal susceptibility, x, is 

The numerator can be evaluated right at q — —Q,u = where it reduces to 
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A(-Q, 0) = (x + _ Q , Q - Xtl- Q ) [(1 - Uxtl Q )(l + Wx^q) + 2U\ + _ p Q qXq P q \ (27) 

Substituting the values for irreducible susceptibilities, we obtain that A = 1/2U and does 
not depend on z. For the total longitudinal susceptibility we then have 



2 



(xi- q r i ^j(q+Q) 2 [i-Y)-2j < 28 ) 

This expression is valid for (tt, Q) phase and for (Q, Q) phase at q x = q y . We see that the 
stiffness for longitudinal fluctuations is pl oc (1 — z/2). For the (ir, Q) phase (1 < z < 2), the 
stiffness is positive, while for the (Q, Q) phase (z > 2) it is negative. This last result implies 
that the homogeneous (Q, Q) phase is in fact unstable. This agrees with our observation in 
the previous section that dE^^ / dx 2 is negative. 

The negative longitudinal stiffness of the (Q, Q) phase was earlier obtained by Dombre 
in the macroscopic calculations in the framework of the Shraiman-Siggia model. He argued 
that this instability leads to a formation of domain walls, but can be prevented by a long- 
range Coulomb interaction. Phase separation at large z is also a possibility [^3|. We however 
have not studied inhomogeneous spin configurations. 

It is also essential to observe that z u behaves as a constant (= z) only at frequencies 
comparable to (q + Q)pf] at larger frequencies, dynamical susceptibility X2d{ui) rapidly 
decreases, and near a spin-wave pole, u 2 ~ 2J 2 (q + Q) 2 , we have X2d(uj) ~ pp\q + Q\/uj <^ 1. 
Then, near the pole, to the lowest order in the hole density we have 

X + - -{^{q + Qf-u 2 )- 1 , (29) 

where the spin- wave velocity, c = \/2J, is the same as in the (tc, tc). This result agrees with 
the macroscopic consideration by Shraiman and Siggia [|I^] and with the Schwinger-boson 



analysis by Gan, Andrei and Coleman 11 . 



B. Transverse spin fluctuations 



We now consider the magnetic susceptibility Xq Z , q associated with the out-of-plane fluc- 
tuations. We found above that this channel is coupled to density and longitudinal spin 
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fluctuations only dynamically. For the full static susceptibility we then have a simple RPA 
formula 

From the above considerations, we expect the Goldstone modes in Xq*-q to be at g = ±Q. 
Consider first the (ir, Q) spiral. Using eq. fl58|) from the appendix A, and expanding near 
q = ±Q, to the lowest order in the hole density we obtain 

where q = n — q and q = n — Q. The second term comes from the integration over the regions 
away from pockets, while the last term comes from the integration inside the hole pockets. 
We see that at q = ±q, i.e., at q — ±Q, \ zz precisely equals 1/2U, and x zz diverges as 
it indeed should. At the same time, neither of the last two terms in (|31|) has a form of a 
quadratic expansion around q = ±Q. To the lowest order in hole density, q = Ux/t, and 
the last two terms in (|5ID cancel each other at any q ~ Q. This means that x zz ls infinite 
in a whole range of momenta which in turn implies that there exists an infinite number of 
other states which are degenerate in energy with the spiral states to the leading order in hole 
density. In appendix B we show that this degeneracy is in fact a quite general phenomenon 
and it exists for an arbitrary ratio of t, t' and U. Furthermore, we found that at least to 
the leading order in t/U, Xq Z - q = 1/2Z7 for q = (7T,q y ) and arbitrary q y . This last result 
means that there exist a whole line of zero modes in x zz - At q = (q x ,q y ), we found that 
Xq Z - q — 1/2Z7 + 0((g x ) 4 ). For typical q x ~ q ~ x, the last term is 0(a; 4 ), i.e., it contains two 
extra powers of x compared to the terms we consider. 

We will discuss the set of degenerate states in the next section, and here merely notice 
that the degeneracy is not related to any kind of broken symmetry and, therefore, should 
be lifted by higher-order terms in the expansion in the hole density which we are proceeding 
to discuss. 

It is not difficult to make sure that the contributions to x zz from the regions far from the 
hole pockets form regular series in powers of q 2 ; odd powers of q disappear due to momentum 
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integration. As typical q ~ q, the next-to-leading order terms have an extra factor of q 2 oc x 2 . 
On the other hand, the expansion near pockets involve only fermions with momenta near 
k = (tt/2, tt/2) and k = (— tt/2, tt/2), i.e., there is no summation over momenta. As a result, 
the next subleading term in x zz from hole pockets has an extra power of x rather than x 2 . 
We computed this term explicitly by expanding in (^) and in the ground state energy (from 
which we extract q) beyond the leading order in hole density, and obtained 



Ux f 7r TO|| + m_i_ 
t \ 8 ^/m\\m± 



(32) 



and 



l- 2 -)-t 



zj q 2 



(33) 



Notice that the correction term in q which explicitly depends on the mass ratio does not 
show up in the expression for \ zz ■ 

We see that the Goldstone mode in \ zz ~ (1 — 2U\ zz )~ l at g = ±q survives as it should, 
but {x zz )~ X stih does not have a form of a quadratic expansion around the Goldstone points 
(Fig. |). For 1 < z < 2, when the (ti, Q) phase is a candidate for the ground state, 1 — 2Ux zz , 
and hence, \ zz '■> is negative in the region q 2 < q 2 . This implies that the (n, Q) planar spiral 
phase is also unstable, and one should look for other candidates for the ground state. 

For completeness, consider also the transverse susceptibility in the (Q, Q) phase. At 
q x 7^ q y , the total transverse susceptibility is positive already at the quadratic order in q: 

l-2U X zz = ^{q x -~q y ) 2 . (34) 

However, along the Brillouin zone diagonal, at q x = q y = q, (x) -1 is again zero, to order 
0(x 2 ), at arbitrary q. We performed calculations to order 0(x 3 ) and obtained 



2U X ZZ = 4x 



3 



H) 


"i _ I _ i£" 




z q 2 



(35) 



The total transverse susceptibility \ zz = X zz I (1 — %Ux zz ) then has one zero associated with 
the Goldstone mode at q — q, and another " accidental" zero at 
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Between these two zeros the total static transverse susceptibility is negative, indicating an 
instability. Contrary to the previously found instability of the (Q, Q) state towards domain- 
wall formation, the latter instability is unlikely to be removed by including the long-range 
component of Coulomb interaction. We however have not performed any further calculations 
for the (Q, Q) phase. 

The result that the planar spiral phase is unstable contradicts the assumption made by 
Shraiman and Siggia that the subleading terms in the expansion over x stabilize the planar 
phase. On the contrary, our results indicate that they do not. 



We conclude this section with a brief consideration of the uniform susceptibility of the 
spiral states. The key point here is that the spiral ordering in the XY plane away from half- 
filling couples the fluctuations along X and Y directions, which at half-filling constituted 



longitudinal fluctuations acquire a correction proportional to the Pauli susceptibility, while 
transverse fluctuations do not. In a spiral phase, these fluctuations are coupled, and as 
a result, the uniform static susceptibility of the ordered state acquires a finite correction 
immediately away from half-filling 



C. Uniform susceptibility 



longitudinal and transverse fluctuations. We have found in Sec. [TI|, that in the (tt, tt) state, 



x + "(o,o) 



1 



X2D 



(37) 



8J 



+ 



4 



where 



X2D = 




rpd _ 

k-Q 



1 



Ei 



k+q 



(38) 
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1 2J 2 

X ™ = 8J T ^ 

We see, therefore, that the step-like correction to the uniform susceptibility is relatively 
small at J <C t. Notice also that fluctuations in the Z direction are decoupled from the XY 
fluctuations, and hence Xoo = (V^-O + 0(x) without any step-like corrections. 

We now proceed to the analysis of the states degenerate with the planar spiral to the 
lowest order in the hole density. 

IV. NON-COPLANAR STATES 

To specify the set of degenerate states, we first return to our results obtained to order 
0(x 2 ) and observe that the zero modes in \ zz are centered around (tt, tt). A zero mode in the 
transverse susceptibility at (tt, tt) implies that the system is indifferent towards generation 
of a spontaneous commensurate antiferromagnetic order along Z direction in addition to 
the incommensurate spin ordering in the XY plane (Fig. [I]). We, therefore, consider a set 
of states having two different SDW order parameters Aj_ = U (S±), and Aii = U(S\\), where 
(S±) and (S\\) are the magnitudes of the off-plane and in-plane components of the on-site 
magnetization, respectively. Performing the mean-field decoupling of the interaction term 
and the diagonalization of the Hubbard Hamiltonian, we obtain 

H = Y,E k (c{c k -4d k ) (40) 

k 

with the energy in the valence (d) and conduction (c) bands given by 



E k = ^(E + + ^A 2 + Elf + Ai (41) 

The ground state energy, to order 0(x 2 ), is given by 

Ux 2 (q*)H 2 , , 

E^ = — -tfx + *±- (42) 

where q* = g(A||/A), and A 2 = A 2 ± + Ajj. For U > t we indeed have A ps U/2. We also 
assumed in © that Ay > x 1/2 . 
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Since we now have two order parameters, there are also two self-consistency conditions. 
The condition on the out-of-plane order parameter Aj_ is: 

£ = ?3B? (43) 

and the condition on the in-plane order parameter Ay is: 

In the limit A^ — > 0, the latter expression reduces to (|TT|) as it should. However, for any 
nonzero An, the two self-consistency conditions have to be satisfied simultaneously, which 
implies that the inverse pitch of the spiral is no longer a free parameter. Specifically, the 
compatibility of the two conditions requires that 

E / E+ = (45) 



However, solving this equation to order 0(x 2 ), we find that 



q* = jx (46) 



which is exactly what one would obtain by simply minimizing the ground state energy fl4"2|) 
with respect to q*. As a result, substituting q* into fl4"2|), we obtain that the ground state 
energy does not depend on Aj_ 

E A± = Ux 2 =E^ (47) 

Clearly then, all states with finite Aj_ are degenerate in energy with the planar spiral, and 
we have to go beyond the leading order in x to see which state actually has the lowest energy. 

The calculations to order 0(x 3 ) proceed in the same way as in the previous section. We 
skip the details and focus on the results. For the inverse pitch of the spiral we obtained from 
the consistency condition (f45f) 



q = — x 



A 2 + A^ 2 



(48) 



Here we again assumed that Ay/A ^> x 1//2 . We see now that the values of q at Aj_ = (eq. 
(j32|) and at A^ — > (eq. fl48|) ) are different to order 0(x 3 ). The reason is, of course, that in 
the case of a planar spiral there is only one self-consistency condition to satisfy. The energy 
of the noncoplanar phase, at An /A ^> x 1 ^ 2 , is given by 



E A± = + ± 



Ux 



t 

where q* is given by (|48|). Let us first discuss the second term in the r.h.s. of (49). This is a 



positive contribution to the energy related to the fact that q* is no longer a free parameter. At 
Aj_ — > 0, the last term in the r.h.s. of (^9|) disappears, and the energy of the noncoplanar 
state turns out to be larger than that of the (tt, Q) state. Clearly then, very close to 
Aj_ = 0, a simple noncoplanar state that we consider is not the best choice. This in fact is 
consistent with the form of the susceptibility in the spiral phase which has a minimum at 
some momentum different from (tt,tt) (see Eq.fl53|)). However, we also see that this energy 
difference is 0(x 4 ), and for all A^ > x 1 ' 2 , the second term in the r.h.s. of fl49D can be 
neglected compared to the third term which is of the order x 3 and negative for z < 2 which 
we consider. This latter term increases with Aj_. Moreover, eq. (|49|) does not contain a 
"restoring force" (i.e., terms ~ x 3 A\). Therefore, by making Aj_ larger and larger, we can 
continuously decrease the ground state energy as long as eq. (f49|) remains valid, i.e., as long 
as Ay/A ^> x 1 / 2 . When A^ nearly reaches A, and Am/A becomes comparable with x 1//2 , 
the expression for the ground state energy of the noncoplanar state becomes more complex. 
In this situation, we found 



Ux 2 Ux 2 



M) 2 + 2(3 2 - ^f(a) 
o 



(50) 



where 



An n q*t A^x 1 / 2 

a = A -I/O . P 



A ± xV2' ^ Ux A,, 
f(a)=(a 2 + ^) -a 3 (51) 



For Ay/A ^> x 1 / 2 , a ^> 1, and expanding in 1/a in (|5lD we recover 
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The solution of the self-consistency condition (f45|) is also more complex for Ay/A ~ x 1 ! 2 . 
We found 



4 

'a 2 H a 



(52) 



At a ^> 1 (but still, Aii <C A), this expression reduces to ([48 



Substituting /3 into (|51~D , we obtain the ground state energy as a function of a single free 
parameter a. The equilibrium value of a (and, hence, of Aj_) can now be obtained by a 
simple minimization of the energy. In a general case, the solution of dE&^/da = is rather 
involved, but for z close to 2, a simple analytical solution is possible. The point is that at 



z ~ 2, the equilibrium value of a is large (~ (2 — z) x / 2 ) 
of 1/a 2 . We then obtain 



E A± = E {w ' Q) + Ux 2 
The energy has a minimum at 



1 / 2 

"a I 1 "" 

cr V 2 



so we can expand (|51|) in powers 



(53) 



24a 4 



72 



' 12(2 -z)' 

For the equilibrium noncoplanar state we thus have 



(54) 



Ea 



Ux 2 
~2z 



(2-z) 



1-^(2-*) 



. / 7 
AW — 



(55) 



12 2-z 

We see that the energy of the noncoplanar state is substantially lower than that of the (it, Q) 
state. What is more, the energy gain in the equilibrium state scales as x 2 , instead of x 3 as 
in (|49|). This implies that the equilibrium state with Aii ~ \fx, strictly speaking, does not 
belong to the original set of degenerate spin configurations, and could be selected already 
in the calculations to order 0(x 2 ). The discovery of the degenerate set of states and of the 
instability of the planar spiral gave us, nevertheless, a hint where to look for the minimum 
of the energy. Notice also that dE,\ ± /dx 2 > 0, i.e., there is no instability towards phase 
separation. 
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Another important point concerns the chirality of the novel state. Although the spin 
configuration in this state is noncoplanar, it is not chiral in the sense that there is no flux 



through a plaquette P^fl . In other words, although, the triple product of three adjacent 
spins along the Y direction, Sij-i ■ (Sij x Sij+i) 7^ 0, the triple product of spins lying in the 
vertices of a minimal triangle is always zero, since all spins along the rows in the X direction 
are parallel to each other. This fact distinguishes our noncoplanar state from the double 
spiral considered in |§b, which has staggered chirality. We emphasize however that, at least 
in the SDW approximation, our noncoplanar state has lower energy than the double spiral. 

V. MAGNETIC SUSCEPTIBILITY OF THE EQUILIBRIUM STATE 

In the presence of a commensurate antiferromagnetism along Z-axis, the equation for 
magnetic instability becomes more complex as now bare susceptibilities with the momentum 
transfer (71*, 7r) are also finite, and the total susceptibility becomes 8x8 problem. In view 
of this, we only computed the susceptibility to the leading order in Aj_. We found that the 
compatibility condition P5| ) of the two self-consistent equations at Aj_ — ^ is equivalent (to 
order 0(x 3 )) to the condition that \ zz diverges at (jr, it), (see Fig H). This extra zero mode 
exists because at Aj_ = 0, the ground state energy as a function of Aj_ has an extremum 
(maximum). We have not performed calculations at Aj_ > Am, but we believe it plausible 
that the spin susceptibility evolve with A^ as it is shown in Fig |^. The equilibrium state 
is an energy minimum (at least, local), and we expect that the static susceptibility of this 
state is positive, diverging only at the Goldstone points. There exists, however, a subtlety 
in determining the locations of zero modes in this configuration, so we will explicitly follow 



the recipe of the Goldstone theorem [^TJ. This theorem states that if J is a generator of 
a symmetry transformation, and the commutator [A, J], where A is some operator, has a 
non-zero average value in the ground state, then the correlator (TA'A) diverges at uo = 0. 
The residue of the quasiparticle pole near the Goldstone point is proportional to [A, J] 2 . In 
our case the corresponding operators and correlators are the following: 
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1. Rotation about the X-axis: J = Sq. If A = S% (tt = (tt,tt))), then [A, J] = —iS, 
Aj_, and therefore Xst'(q) diverges at q = 7T. The residue of the pole is proportional to 
Aj_. If instead we choose A = S±q, then [A, J] = iS±Q, and we find a divergence in 
Xg t z (q) at q = ±Q; the residue of the pole is proportional to Ajj. 

2. Rotation about the y-axis: J = Sq. This case is analogous to the previous one. The 
divergences occur in Xst( n ) with the residue of the pole oc A\, and in Xg t z (±Q) with 
the residue oc Ajj. 

3. Rotation about the Z-axis: J = Sq. Choosing A = Sq and A = SZq we find diver- 
gences in Xst~(Q) an d Xst + ( — Q)i i n both cases the residue of the pole is proportional 
to Ajj. 

Combining these results, we find that in the equilibrium noncoplanar state, the in-plane 
static susceptibility x + ~ m fact has two poles, q = (tt, 7r) and the other at q = —Q. 

Then, for q not too far from (tt, tt), this susceptibility can be approximated as 



^ W ?^ + ^> (56) 
(q-vr) 2 (q+Q) 2 



where the residues x-k and Xq are proportional to A^ and Ajj respectively. Because Ay ~ t/x, 
the residue of the pole at the incommensurate wave vector q = — Q is suppressed with respect 
to the pole at the commensurate wave vector (tt, tt), and the form of the static susceptibility 
is very similar to that for the (tt, tt) state (Fig. 

The out-of-plane static susceptibility x zz also has two poles 

(q-Qr (q + Q) 

but here both poles are suppressed as the residue, Xz, is proportional to Ajj, and, therefore, 
to x. Again, the form of the susceptibility is very similar to that for the (tt, tt) state. On the 
contrary, in the planar spiral state, both transverse and longitudinal susceptibilities have 
Goldstone modes at q — ±Q with the residue of the pole proportional to the total on-site 
magnetization (Fig. |). 
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VI. CONCLUSIONS 



Here we summarize the main results of this paper. We used spin-density-wave formal- 
ism and studied various magnetic phases of the 2D Hubbard model at low doping when 
a long-range magnetic order is still present. We found that the equilibrium spin configu- 
ration depends on the value of the dimensionless parameter z = 4Tx 2D , where T is the 
scattering amplitude of two holes (T — U in the mean-field approximation), and x 2D is 
the Pauli susceptibility of holes which at low doping occupy pockets located at (7r/2,7r/2) 
and symmetry related points in the Brillouin zone. In 2D, Pauli susceptibility does not 
depend on the carrier concentration: x 2D = y /fn\\m_ s _/2TT. We found that for z < 1, the 
commensurate antiferromagnetic state is stable. For z > 2, the spiral (Q, Q) phase is the 
equilibrium configuration at the mean-field level, but we found that this configuration has 
negative longitudinal stiffness and therefore is unstable against domain wall formation. This 
result agrees with the macroscopic analysis in ||. The intermediate case, 1 < z < 2, is the 
most interesting from a theoretical point of view. Here we found that the equilibrium con- 
figuration at the mean-field level is a (tt, Q) spiral introduced by Shraiman and Siggia. The 
longitudinal stiffness in this configuration is positive, but the transverse stiffness vanishes 
to the leading order in hole density. This in turn implies that the spiral phase is degenerate 
in energy with many other spin configurations, and the equilibrium state only appears as 
an "order from disorder" effect. We performed calculations beyond the leading order in 
the hole density and found that the equilibrium state is not a planar spiral but rather a 
noncoplanar spin configuration which contains both, (it, tt) antiferromagnetism along one 
direction in spin space, and (tt, Q) spiral in the orthogonal plane. The latter result suggests 
a novel scenario of spin reorientation with doping for 1 < z < 2, different from the one 
suggested by Shraiman and Siggia. In their picture, upon doping spins remain in the same 
plane as at half-filling, but twist into a spiral with incommensurate momentum (tt,Q). In 
our scenario, the commensurate antiferromagnetic ordering (same as at half-filling) does 
not vanish, as doping only introduces a transverse component of the order parameter which 
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forms a spiral in the plane perpendicular to the direction of the commensurate order. This 
transverse component is small to the extent of x, and the low-T behavior at finite doping 
remains nearly the same as in the commensurate antiferromagnet [^,Q(Fig. |2[). 

It is essential that our analysis has been performed only for frequencies smaller than the 
energy scale AE associated with the lifting of the degeneracy. At larger frequencies, the 
static selection may be irrelevant, and one has to solve the full dynamical problem which 
presents a technical challenge. 

The above analysis is valid for the magnetically ordered phase. We therefore cannot 
pretend to resolve the known discrepancy between neutron scattering and NMR experiments 



in La 2 - x Sr x CuO^ [p7| , p8|| , both of which have been performed well inside the metallic phase. 
We merely note that, as neutron data indicate, the incommensurability at (it, Q) observed 
at 7.5% and 14% doping is not correlated with the magnetic behavior in the ordered phase. 
This implies that our result that in the ordered state, susceptibility is always peaked at 
(7r,7r), does not contradict the neutron data. 

It is our pleasure to thank A. Abanov, V. Barzykin, P. Chandra, P. Coleman, R. Gooding, 
R. Joynt, D. Khomskii, D. Pines, S. Sachdev, Q. Si and A. Sokol for useful discussions. The 
work was supported by the University of Wisconsin-Madison Graduate School and Electric 
Power Research Institute. 



VII. APPENDIX A 

In this appendix we present the results for the irreducible susceptibilities in the spiral 
phases. Each of the susceptibilities below was obtained by the standard SDW manipulations. 
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Further, 
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Next, 
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Again, upper sign is for x +z , lower for Finally, 
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(63) 



We also have - X q + 2 Q,- q ~ X q * {q+2 Qy xJ+o^( 9 +2Q) - X g +2<f-( g +Q)> xJ-J+Q) ~~ X 9+S,-? 
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VIII. APPENDIX B 



In this appendix, we show that the zero modes in the transverse susceptibilities of the 
two spiral phases exist, to order x 2 , independent of the ratio oit/U and t'/U. To see this, 
consider the transverse susceptibility right at q — (ir,Tr). We will show that \ zz = 1/2Z7, 
i.e., the total X^ z —k = X zz (^ ~ 2?7x z ^) _1 diverges despite the fact that for the spiral states, 
(71", 7r) is not the ordering momentum. For definiteness, we will perform the calculations for 
the (71", Q) phase. The calculations for the (Q, Q) phase proceed in the same way, and the 
final result is valid for both spiral states. 

Expanding in ([58]) to the second order in q ~ x we obtain 

111 A 2 1 

zz \ ^ \ ^ f — — \ 2 

1 A 2 1 

8iV E fJ ek + ^ (MY Ei - Ei +n (64) 

fc + 7T 

where 7r should be interpreted as 2D momentum (7r, 7r), and as before, q = tt — Q. It is also 

convenient to redefine the momenta such that e^ - = {tk+Q i e fc)/2. For the (vr, Q) spiral, 
we then have 

6fe + efc +vr = 4|t'| cosfc x sinfc^ g, (65) 

and also 

4+n ~4 = ~ 2t V sink V ( 66 ) 



Substituting (^) and flSBD into (0), we find after some simple algebra 

1 (sin A:^) 2 _ {t') 2 f ^ (cosfc x sin k y ) 2 _ tqx_ 

Xn '- n 2U 4iV Y (£ fc ~) 3 AT Y ( E fc") 3 452 ' 

In obtaining this result, we used the relation 

1 _ / 1 1 \ 9/^2=2 1 \- {sinky cosk x ) 2 (e^) 2 
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which can be derived by straightforward computations using fl65|) . In (^7]) and (|68|), the 
summation is over the whole Brillouin zone. 



Notice that the pocket contribution (a term ~ tqx in (|6?|)) is the same as in the analysis 
in the bulk of the paper, where we assumed that t', t <C U . This is simply related to the fact 
that the pockets are located at (vr/2, ±7r/2) where both e^ande^ are small compared to A 
independent of the ratio of the parameters. 

We now need the exact relation between q and x, valid to the first order in x, but for 
arbitrary t, t' and U. To find this relation, we again compute the ground state energy of 
the (it, Q) spiral, but this time without assuming that U is large compared to the hopping 
integrals. Doing the same computations as in Sec. |J, we find 

£j(7r,Q) = _ t y X _|_ ^2-2 2 cos 2 k y - sin 2 k y ^ sin 2 k y (e k ^ 2 1 



2(0 2 A¥^E Sin2 ( ^f fcx - (69) 
The equilibrium q then satisfies 
tq X _ ^2-2 J_ ( 2 cos 2 k y — sin 2 k y ^ sin 2 k y {e k ) 2 \ _ 4.^f;') 2 A 2 q 2 — ^ S * n2 ^ y C ° s2 ^ x 



Substituting this result into (0), we obtain 



(70) 



= 2!7 + 5 tV A ' (71) 



where 



a 1 A 2 sin 2 k y 1 cos 2 fc y 

Notice that all terms with t' are cancelled out. Finally, to evaluate \ zz to order x 2 , we 
actually need A only for q = 0. In this case, = — 2t(cosk x + cosfcy), and integrating by 
parts in (|72|) , we immediately obtain that A = is independent of the ratio of t/U. We thus 
find 
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xlU=yj + 0{x*) (73) 

This result implies that the zero modes in the transverse susceptibility exist at an arbitrary 
ratio of the parameters of the Hubbard model provided that the magnetic ordering at half- 
filling is commensurate, and doped holes form pockets at (n/2,n/2) and symmetry related 
points. 
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FIGURES 

FIG. 1. Spin configuration of a noncoplanar state. Arrows with thick ends point out of the 
plane, while those with thick tails — into the plane. This configuration is different from the double 
spiral considered in |^]b. 

FIG. 2. Two adjacent spins in the equilibrium configuration. The in-plane component, 
S± ~ x 1//2 , is small compared to the off-plane component, S\\. 

FIG. 3. The out-of-plane static susceptibility for the planar (tt, Q) spiral. The susceptibility is 
negative around (tt, it) indicating instability towards spontaneous magnetization in the out-of-plane 
direction 

FIG. 4. The in-plane static susceptibility for the equilibrium noncoplanar state, and a non- 
coplanar state with vanishing Aj_. 

FIG. 5. The out-of-plane static susceptibility for three noncoplanar states: one with A^ — ► 0, 
another with an intermediate value of Aj_, and the third with the equilibrium value of Aj_. 
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t Oc+-(q;G) = 0))- 1 



equilibrium state 




